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ABSTRACT 

We present optical and near-infrared photometry, as well as ground-based 
optical spectra and Hubble Space Telescope ultraviolet spectra, of the Type la 
supernova (SN) 2001ay. At maximum light the Si 11 and Mg 11 lines indicated 
expansion velocities of 14,000 km s~^, while Si III and S II showed velocities of 
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9,000 km s~^. There is also evidence for some unburned carbon at 12,000 km s~^. 
SN 2001ay exhibited a dechne-rate parameter Ami5(i?) = 0.68 ± 0.05 mag; this 
and the 5-band photometry at t > +25 d past maximum make it the most slowly 
declining Type la SN yet discovered. Three of four super-Chandrasekhar-mass 
candidates have decline rates almost as slow as this. After correction for Galactic 
and host-galaxy extinction, SN 2001ay had Mb = —19.19 and My = —19.17 
mag at maximum light; thus, it was not overluminous in optical bands. In near- 
infrared bands it was overluminous only at the 2a level at most. For a rise time 
of 18 d (explosion to bolometric maximum) the implied ^^Ni yield was (0.58 ± 
0.15)/a Mq, with a = L^s,^/ E^i probably in the range 1.0 to 1.2. The ^^Ni yield 
is comparable to that of many Type la supernovae. The "normal" ^^Ni yield and 
the typical peak optical brightness suggest that the very broad optical light curve 
is explained by the trapping of the 7 rays in the inner regions. 

Subject headings: supernovae: general — supernovae: individual (SN 200 lay) 



Introduction 



Phillips! (Il993[ ) first established that Type la supernovae (SNe) are standardizahle candles 
at optical wavelengths: there is a correlation between their absolute magnitudes at maximum 
light and the rate at which these objects fade. This fact allowed Type l a SNe to be used 



to determine that the expansion of the Universe is currently accelerating (IRiess et al.lll998 



Perlmutter et al.lll999l ). More recently, we hav e discovered that in near-infrared (IR) photo 



metric bands Type la SNe are standard candles (Meikle 200o[ 



2004al : IWood-Vasev. et al. 1 120081 iFolateUi et al. 



2010 



Krisciunas. Phillips, fc Suntzefl 



Mandel. Narayan. fc Kirshnerll201ll ). 



Except for a fraction o f the rapidly declining T ype la SNe whose prototype is SN 1991bg 
( iFilippenko et al.lll992t iLeibundgut et al.lll993l ). which peak in the near-IR after the time 
of S-band maximum, the near-IR absolute magnitudes at maxin aum light are at most 



only slightly depen dent on the decline-rate parameter Ami5(-B) (IKrisciunas et al 
Folatelh et ZIboiol ). 



2009 



We do not understand well how to model a Type la SN. For much of the previous 
decade we thought that (most) Type la SNe were carbon-oxygen white dwarfs that approach 
the Chandrasekh ar limit (1.4 Mq) owing to mass transfer from a nearby companion (for a 



re view, see 



_e^ 



Livio I I2OOOI ). Some Type la SNe might be mergers of two white dwarfs 



( llben &: Tutukovlll984l ). Over the past few years the double-degenerate scenario has gained 
prominence among many researchers. Unfortunately, it is not yet possible to determine 
with a high degree of probability that a given Type la SN is a single or double-degenerate 
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event, except perhaps for the few that produced more than a Chandrasekhar mass of fusion 
products; these would result from double-degenerate mergers. An excellen t sura mary of 
evidence for different kinds of progenitors of Type la SNe is given by iHowelll (120 111 ). 



The spectroscopic classification scheme stipulates that Type I SNe do not ex hibit obvious 



hydro gen in their spectra, an d that Type II SNe do show prominent hydrogen ([Minkowski 



194ll ): see iFilippenko I ( 1l997l ) for a review. Type la SNe show Si II in absorption, blueshifted 
from its rest wavelength at 6355 A. Two objects have been s hown to have str ong Si II in ab 



sorption and Ha in emiss ion: SN 2002ic (IHamuy et al.ll2003l ) and SN 2005gj (lAldering et al. 



20061 : iPrieto et al.ll2006a( ). The interpretation is that these objects are Type la SNe interact- 
ing with circumstellar material (CSM) or (less likely) with the general interstellar medium 
(ISM). 



SN 2006X was found to exhibit variable Na I D absorption (IPatat et al.l 12(3071). Two 
other exam ples showing variable Na I absorption are SN 1999cl (IBlondin et al.l l2009l ) and 
SN 2006dd (IStritzinger et al.ll2010f l. Light from the SNe is ionizing the circumstellar medium. 
This leads to recombination and line variations. It should be pointed out that SNe 1999cl 
and 2006X have large reddening that m ay involve multiple scattering as well as normal dust 
extinction dWang lEoOsI : ICoobailboOsl ). 



It is possible to relate the peak bolometric luminosity of a Type la SN and the implied 
amount of radio active ^^Ni produced in the explosion. Most events generate between 0.4 and 
0.7 Mq of ^^Ni Jstritzinger et al1l2006h . The remarkable SN 2003fg (also known as SNLS- 
03D3bb; Howell et al. 20061 ) was sufficiently overluminous that it was regarded to have been 
caused by the explosion of more than 1.4 M© of carbon and oxygen. Th ree additional super - 



Chandrasek! 



lar-mass candida t es have recently be en found: SN 2006 gz (IHicken et al.l 



SN 2007if (lYuan et al. 



Tanaka et a. 



Scalzo et al 



2010 



2010; 



Silverman et al. 



Scalzo et al.l 2010), and SN 2009dc (Yamanaka et al. 



2010[ iTaubenberger et al.l l201l[ ) As summarized by 



2007) 



2009 



2010l . Table 4), these objects can be up to a magnitude more luminous at optical 
wavelengths than typical Type la SNe. Furthermore, three of these four have extremely 
slowly declining light curves (Amis (5) ~ 0.7 mag). In a double-degenerate system, the 
more massive white dwarf can tidally disrupt the less massive white dwarf, which has a larger 
radius. A disk or shell of material is created, retarding the expansion of the more massive 
white dwarf once it explodes. Another objec t with a very slow decline rate was SN 200 5eq, 
a target of the Carnegie Supernova Project ( iFolatelli et al.ll2010l : IContreras et al.ll2010l ). 



In this paper we present optical and near-IR photometry, plus a plethora of spectra, of 
the equally remarkab le SN 2001ay. Preli minary analysis of the light curve and spectra of 
SN 2001ay is given by lPhillips et al.l ( 120031 ). Two Hubble Space Tele scope (HST) STIS spectra 
(from 2001 May 2 and May 9) have already been published by iFoley. Filippenko. fc Jha 
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( 120081 ) ■ in a study of possible luminosity indicators in the ultraviolet (UV) spectra of Type 



la SNe. 



2. Optical and IR Photometry 



SN 2001ay was d iscovered on 2001 April 18.4 (UT dates are used throughout this paper) 
bv ISwift fc LjfcoOlh . Its position was a = 14^'26°>17.0^ 5 = +26°14'55.8" (J2000), some 
10.3" west and 9.3" north of the nucleus of the spiral galaxy IC 4423. Basic information 
on SN 2001ay is given in Table [TJ Figure [T] shows the host galaxy and SN, and identifies a 
number of field stars of interest. 

A sizable fraction of our optical photometry of SN 2001ay was obtained with the 0.76- 



m Katzman Automat ic Imaging Telescope (KAIT) at Lick Observatory (ILi et al.l 12000 



Filippenko et al.ll200ll ). KAIT images are 6.7' x 6.7' in size, with a scale of 0.80" per pixel. 
For the first three nights of photometry with the Lick Observatory Nickel 1-m telescope, the 
CCD camera had a chip with 0.30" pixels in 2 x 2 readout mode and gave a field size of 
5.1' X 5.1'. Thereafter, a new chip was installed, having 0.37" pixels in 2 x 2 readout mode 
and a field size of 6.3' x 6.3'. The two different CCD chips employed with the Nickel telescope 
had considerably different quantum efficiencies at blue wavelengths. 

The UBVRI photometry of some field stars near SN 2001ay is found in Table [2l Cali- 
bration of optical photometry was accomplished from observations carried out on five pho- 
tometric nights (two with the CTIO 0.9-m, two with the N ickel, and one with the CTIO 
1.5-m telescopes) using the standard stars of Landolt (Il992). A subset of o ur photometry 
(the KAIT data) has already been published by iGaneshalingam et al.l (120101 . Table 6), but 
without K- and S-corrections. Our data presented here were redu ced in the IRAF environ- 
mentl^ Transformation of the d ata to the sy s tem of Landolt (1992) was accomplished using 
equations equivalent to those of IWang et al. J2009ah . In their Table 1, the reader will find 
characteristic color terms used for the transformations. Our own direct determinations of 
these color terms are entirely consistent with their values. 

Table [3] gives the Js and if-band photometry of some of the field stars. We give the 
averages obtained on seven photometric nights using the Las Campanas 1-m Swope telescope. 
The calibration was accomplished using a stand-alone version of daophot, some FORTRAN 



^^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation (NSF). 
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programs written by one of us (N.B.S.), and the IR standards of iPersson et al.l (Il998[ ). Thus, 
our near-IR photometry of SN 2001ay is on the iPersson et al.l (119981 ) photometric system. 



We can compare the photometry of the principal I R secondary standard ( "Star 1" ) with 
data from the Two-Micron All Sky Survey (2MASS; ISkrutskie et al.l 120061 ). Whereas we 
obtained J, = 13.718 ± 0.014 and H = 13.251 ± 0.011 mag, 2MASS found J = 13.718 ± 
0.027 and H = 13.214 ± 0.030 mag. These values are statistically in agreement. For "Star 
6" we obtained J, = 15.153 ± 0.021 and H = 14.881 ± 0.008 mag, whereas 2MASS found 
J = 15.132 ± 0.045 and H = 14.807 ± 0.060 mag. 

Tables m and [7] give the optical and IR photometry of SN 2001ay, respectively. Table H] 
lists point-spread function (PSF) magnitudes from KAIT, the first three Nickel 1-m nights, 
and one night of the CTIO 1.5-m telescope; the subtracted template images were obtained 
with KAIT on 2002 July 9 and 10. For the single night of imaging with the CTIO 1.5-m 
telescope, the ?7-band magnitude was determined from aperture photometry; no t/-band 
template was available. For the final night of Nickel 1-m photometry, PSF magnitudes 
were derived, but without template subtraction. Aperture photometry without template 
subtraction was conducted for all other optical images|lf| 



To compare the photometry of different SNe, one shou 
frame equivalent by subtracting the so-called K-corrections ( 


d transform the data to the rest- 


Oke & Sandaeelll968: 


Hamuv e 


t al. 


1993; 


Kim. Goobar. & Perlmutter 


1996 


Horn et al. 


2002 


Nueent. Kim. & Perlmuttei 


2002 


)• 



For each SN one computes the time since maximum brightness, and then obtains a rest-frame 
timescale by dividing the differential times by 1 + 2;, where z is the object's heliocentric red- 
shift. We h ave wor k ed ou t the optical K-corrections using actual spectra of SN 2001ay, 



adopting the iBesselll (Il990l ) filter profiles as reference; these are given in Table [5] and shown 
in Figure m 

To account for differences in SN photometry that result from the use of different 
telescopes, CCD chips, an d filters, we use t he na ethod of spectroscopic corrections (S- 
corrections), as outlined by IStritzinger et al.l (120021 ) and iKrisciunas et al.l (120031 ). For op- 
tical bands we adopt as reference the filter specifications of IBesselll (Il990l ). At minimum 
one constructs the effective filter profiles using the laboratory transmission curves of the 
filters multiplied by an appropriate atmospheric function that includes atmospheric absorp- 
tion lines; then, one multiplies that result by the quantum efficiency of the CCD chip as a 
function of wavelength. Telescope and instrument optics also contribute to the effective filter 



^^Due to the long delay in assembling the data and writing this paper, some images are no longer retrievable 
from magnetic storage media. Only aperture photometry was possible in 2001, when some of the data were 
reduced, owing to the lack (at the time) of host-galaxy templates. 
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profiles. For the KAIT and Nickel 1-m corrections we have used an atmospheric functioio 
appropriate to Lick Observatory (elevation 1290 m). For the CTIO 0.9-m data we have used 
a different atmospheric function appropriate to Cerro Tololo (elevation 2215 m). In practice, 
one arbitrarily shifts the effective filter profiles toward longer or shorter wavelengths so that 
synthetic photometry based on spectra of standard stars gives color terms that match those 
obtained from actual photometry of ph otometric standards. We used spectra of 50 stars from 



the sample of IStritzinger et al.l (120051 ) to calculate our synthetic magnitudes using an IRAF 
script written by one of us (N.B.S.). This script was then used to calculate the S-corrections 
based on spectra of SN 200 lay itself. 

In Figure [3] we show the B- and V-band S-corrections for SN 2001ay. One adds these 
corrections to the photometry. As one can see in Figure [31 if an object like SN 2001ay 
were observed with the CTIO 0.9-m and Nickel 1-m telescopes, there could be photometric 
discrepancies as large as ~ 0.06 mag in the B band, depending on the epoch. In the V band 
the corresponding differences are much smaller. 

Even after the application of the S-corrections, we find that the Nickel 1-m 5-band 
data based on aperture magnitudes are systematically 0.09 mag fainter than the KAIT 
measurements. For the first three nights of Nickel imaging we were able to eliminate these 
systematic differences by means of template subtraction. The images subsequently obtained 
with the Nickel 1-m telescope and its newer chip did not lend themselves to image subtraction 
using the KAIT templates. (Depending on the number of field stars, the angular size of a 
nearby galaxy in the field, and the seeing, the image rescaling and remapping algorithm can 
fail.) As a result, we have devised a third set of corrections to the photometry based on images 
with the new chip in the Nickel camera. Using addstar within the daophot package of 
IRAF, we were able to add an artificial star to the KAIT image-subtraction templates at the 
pixel location of the SN. This artificial star can be scaled to give a standard magnitude of a 
desired value such as B = 17.97 and V = 17.14 mag, just about the brightness of the SN on 
2001 May 14 (JD 2,452,043.8), in the middle of the run of Nickel 1-m aperture photometry. 
Artificial stars having the identical PSF magnitude as the fake star at the location of the 
SN are placed at blank places in the KAIT templates using addstar. Since the Nickel 1-m 
aperture photometry was typically derived with a software aperture of radius 10 pixels (px) 
and a sky annulus ranging from 12 to 20 px, we can then compare the aperture magnitudes 
of the fake SN with aperture photometry of the fake stars in the blank places of the image. 
We use a software aperture radius and sky annulus of identical size in arc seconds to that 
used for the Nickel 1-m photometry. In this way, we can obtain an estimate of the systematic 



^^The atmospheric functions combine the mean broad-band extinction values with atmospheric absorption 
lines. 
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errors of the SN aperture photometry. 

We determined that the i?-band magnitudes from Nickel 1-m aperture photometry of 
the SN were 0.08 mag fainter than what we would have measured without the presence of 
the host galaxy. Similarly, in the V band we found that the Nickel 1-m aperture magnitudes 
were fainter by 0.05 mag. For the two nights of CTIO 0.9-m aperture photometry, when the 
SN was more than 1 mag brighter, the SN was too faint by 0.03 mag in B and V without 
this correction. Thus, by applying the S-corrections and these additional corrections, we 
can reconcile almost all systematic differences in the B- and l^-band photometry obtained 
with different telescopes, using different cameras, and using different data-reduction methods 
(i.e., PSF magnitudes with image-subtraction templates vs. aperture magnitudes without 
subtractions). 

Further justification of this method comes from a consideration of the optical photome- 
try from the first Nickel 1-m chip and the photometry from one night with the CTIO 1.5-m 
telescope. Data from these four nights can be derived using the KAIT host-galaxy templates. 
In the B band, photometry of the SN using PSF magnitudes and image subtraction was 0.07 
to 0.16 mag brighter than aperture photometry using an annulus for the subtraction of the 
sky level. In the V band, the photometric values of the SN on the first three Nickel 1-m 
nights were 0.06 mag brighter than values derived from aperture photometry. This result 
is contrary to typical experience; normally, aperture photometry of a SN is brighter than 
expected (not fainter) if the SN is located on top of some part of the host galaxy, because 
the light in the aperture is not entirely due to the SN. That this is not the case here must 
be due to the relative distributions of host-galaxy light at the SN position and in the sky 
annulus. 

Differences between aperture photometry and PSF photometry are usually larger in 
redder bands. In fact, similar experiments with adding artificial stars to the KAIT templates 
of SN 2001ay show that the aperture magnitudes with the Nickel 1-m telescope are too bright 
by 0.01 mag in R and I for the size of the aperture and annulus used. We never obtained 
near-IR subtraction templates with the camera used on the Las Campanas Observatory 1-m 
telescope, and that camera has since been decommissioned. So, we must adopt the available 
near-IR aperture magnitudes. 

The J-band photometry is particularly problematic from 15 to 30 d after the time of 
-B-band maximum. Consider the effective filter profiles shown in Figure HI While the KAIT 
J-band filter very well approximates the Bessell / filter, the filter used on the Nickel 1-m 
telescope does not. As the SN achieved its reddest optical colors, a significant excess amount 
of light is let through by the Nickel 1-m /-band filter. This leads to positive S-corrections; 
to correct such photometry to Bessell filter photometry requires making the Nickel 1-m 
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photometry fainter. Our attempts to reconcile up to 0.35 mag differences between the 
KAIT and Nickel 1-m photometry were not successful. Apparently, the effective filter profile 
of the Nickel 1-m camera is even more nonstandard than our profile based on laboratory 
traces of the filter, knowledge of the quantum efficiency as a function of wavelength, and a 
generic atmospheric function applicable to Lick Observatory. While we list all of our optical 
photometry in Table HI the Nickel 1-m J-band photometry is not included in our plots or 
used in the analysis. 

The R- and J-band S-corrections are shown in Figure O Note that the photometry in 
Table S] includes the K- and S-corrections, plus the corrections mentioned above for CTIO 
0.9-m and Nickel 1-m photometry reduced using aperture magnitudes. Since the K-corrected 
photometry from one night with the YALO 1-m telescope and four nights with the LCO 2.5- 
m telescope are in good agreement with photometry obtained with other telescopes, and 
we have no effective transmission curves for the LCO 2.5-m filters used, we have derived 
no further corrections for this small fraction of our photometry. The interpolated BVRI 
K-corrections and S-corrections are listed in Table [6l^i 

In spite of the rationale outlined above to reconcile as much of the optical photometry 
as we can, for a derivation of the maximum brightness and decline rate of SN 200 lay we 
restrict ourselves to the photometry based on PSF magnitudes with image subtraction. 
After subtracting the derived K-corrections and adding the S-corrections, we scale the time 
since maximum light to the rest frame by dividing by 1 + z. We derive a time of S-band 
maximum of JD 2,452,022.49 (= 2001 April 23.0), with an uncertainty of perhaps ± 0.8 d. 
The decline-rate parameter found from our best S-band data is Ami5(-B)= 0.68 ± .05 mag. 



For comparison SNe 2 005eq had Ami5(E) = 0.72 ± 0.02 mag (IFolatelh et al.ll2010f l. though 



Contreras et all (l2010h give A mi JEW 0.78 ± 0.01 mag. SN 2009dc had Ami5(5)= 0.72 



± 0.03 mag (ISilverman et al.ll2010f ). SNe 2001ay, 2005eq, 2006gz, 2007if, and 2009dc had 
extremely slow decline rates. Since it is well established that the energy budget of a Type 
la SN is related to the decline rate, we naturally wonder if SNe 2001ay and 2005eq are also 
super-Chandrasekhar-mass candidates. 

The B- and \^-band light curves of SNe 2001ay, 2005eq, and 2009dc are illustrated 
in Figure O For the first three weeks after maximum light SNe 2001ay and 2009dc were 
remarkably similar, but at t ~ 25-40 d the light curves diverge. As the peak-to-tail ratio in 
the B- and V^-band light cu rves sheds light on the central density and progenitor mass of 



single-degenerate explosions (IHoflich et al.ll2010l ). a comparison of the observational features 



^"^At present, pipelines for SN surveys use template spectra at different epochs since maximum light, warp 
them, and, adopting appropriate effective filter profiles, calculate the K- and S-corrections at the same time. 
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and modeling of very slow decliners may yield similar insights. 



In Figure 
locus of Lira 



7| we show the unredden e d B — V colors of SN 2001ay, the unreddened 
19951 ) and [phimps et al.l f ll999h . and the Amis (5) = 0.83 mag locus of 



Branch et al. 



Prieto. Rest, fc Suntzefj J2006bh . On the basis of some of the spectra published here, 



(120061) c l assified SN 2001ay as a "broad-line" Type la SN. This is the sub- 



type that IWang et al.l (l2009bl ) suggest is intrinsically redder than normal Type la SNe, or 
occurs in dustier environments, with Ry ^ 1.7. However, Figure [7] shows that, if anything, 
SN 2001ay is bluer than other normal Type la SNe. 

In Figure [8] we show the corresponding R- and /-band light curves of SN 2001ay, along 
with loci derived from the data of SNe 2005eq and 2009dc. The behavior of SN 2001ay in 
these photometric bands is clearly more like that of SN 2009dc than of SN 2005eq. Note that 
the /-band secondary maximum of SN 2001ay is essentially as bright as the first maximum. 

Figure [9] shows the near-IR photometry of SN 2 00 lay and t he other slow decliners , 



SNe 2005eq (IFolateUi et al. 



SN 1999aw fIStroker et al 



■■20101 : IContreras et al.ll2010f ) and 2009dc fIStritzinger et al.ll201lD . 
2002h and SN 2001ay were the first objects known to exhibit such 



flat //-band light curves. (Both //-band light curves are admittedly somewhat ragged, how- 
ever.) SN 2005eq was the first to show such early near-IR maxima {t < —7 d). Interestingly, 
the //-band brightness of SN 2009dc increased steadily over the time frame that the bright- 
ness of SN 2001ay was constant. 



3. Spectroscopy 

We obtained spectra of SN 2001ay on 20 dates using 9 different telescopes; 7 of those 
nights had spectra with more than one telescope. Spectra were obtained with the 1.5-m 
telescope of the Fred L. Whipple Observatory (FLWO) on 11 nights. Eight spectra were 
taken on four nights with the MMT 6.5-m telescope; four with the 2.16-m telescope at 
the Xinglong Station of the Beijing Astronomical Observatory (BAO); three with the Las 
Campanas Observatory 2.5-m du Pont telescope; two each with the Lick 3-m Shane reflector 
(Kast spectrograph) and HST; and one each with the Kitt Peak 2.1-m, Kitt Peak 4-m, and 
Keck II telescopes. The spectroscopic observations with the FLWO 1.5-m and the MMT are 
summarized in Tabled we may refer to this as the "CfA set." A log of other spectroscopic 
observations is given in Table M Some, bu t not all, of the data were obtained at (or close 



to) the parallactic angle (iFilippenkd Il982l ) in order to minimize the effects of atmospheric 
dispersion. 



In Figure [12] we show a temporal sequence of spectra. Here we have combined the May 2 
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spectrum from HST with the April 29 spectrum from the KPNO 2.1-m telescope. The rather 
noisy BAO spectra from May 20, 25, and 30 are not illustrated, but the May 10 spectrum 
from BAO fills a gap in the temporal coverage. We show additional spectra in Figure [T31 



Figure [H] includes a portion of our spectrum of SN 2 001ay obtained with the Keck 
Echellette Spectrograph and Imager (ESI: ISheinis et al.ll2002l ) from 2001 April 22 {t = —1 d). 
We find that this maximum- light spectrum exhibits S II and Si III at —9000 km s~^. The 
Si II line (rest wavelength Aq = 6355 A) and the Mg II line (Aq = 4481 A) are blueshifted 
by 14,000 km s~^ at this epoch. There is some gas of these species blueshifted as much as 
20,000 km s~^. As shown in Figure [T2| hj t = +17 d the Si II absorption is only blueshifted 
about 10,000 km s'^H 

In our Keck spectrum there is a small absorption dip at an observed wavelength of 
A ~ 6520 A. If this is due to gas in or near the SN or gas in the ISM of the host galaxy, 
the wavelength of this line in the frame of the host galaxy is 6315 A. This may be due to 
C II (Ao = 6580 A) blueshifted by 12,000 km s~^. However, other explanations are possible. 
An absorption feature at 6520 A in our Galaxy or in the Earth's atmosphere could lead 
to a misidentification, given that telluric lines were not removed from the Keck spectrum. 
One possible such feat ure is telluric absorption at 6515 A due to atmospheric water vapor 
( iMatheson et al.ll2000l . Appendix). We have better evidence of C II absorption in SN 2001ay 
from the Aq = 4745 A line, which is seen in the galaxy's rest frame at 4555 A, corresponding 
to an identical blueshift of 12,000 km s~^. We see no evidence of the corresponding C II line 
with Ao = 7234 A. 

One might expect the C II to be concentrated in the outer layers, and therefore at 
a higher velocity o f approach than sulfur or silicon (IKhokhlov. Miiller. Sz Hofiichl Il993l : 
Gamezo et al.l 120031 ). But if there is mixing or dumpiness, most species could be observed 
over a broad range of velocities. 



In Figure [15] we see the two components of Na I D from gas in the Milky Way and 
gas in the host galaxy of SN 2001ay. The presence of these lines implies some nonzero 
amount of reddening and extinction. We find equivalent widths of 0.090 ± 0.018 A for the 
Galactic sodium l i nes an d 0.210 ± 0.018 A for the host-galaxy sodium lines. Using the 
Munari fc Zwitterl ( 119971 ) calibration of the equivalent width of the sodium lines with B — V 
color exce ss, we find E(B - V)a,^ = 0.026 ± 0.006 mag and E{B - 1/)host = 0.072 ± 0.008 
mag. The lMunari fc Zwitten (119971 ) calibration shows a scatter of ±0.05 mag for E{B — V), 



which is a more realistic estimate of the uncertainty of the reddening toward SN 200 lay. 



For the rest wavelengths of these and other species, see Table 1 of I Wang et all toodl) . In their table, 
however, the rest wavelength of Mg II should read 448.1 nm, not 447.1 nm. 
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^ The Galactic component may b e compared to E{B 
Schlegel. Finkbeiner. fc David ( 1l998l ). 



V) = 0.019 mag obtained by 



In Figure [TBI we see the (blueshifted) velocities of Si II A4130 and A6355, with respective 
error bars of ±420 and ±330 km s~^. The A6355 line exhibits the larger velocities because its 
opacity is greater; we are measuring material farther out in the expanding fireball. The high 
velocities suggest that the outflow of SN 2001ay is essentially unimpeded, which i s consistent 
with it having been a single-degenerate explosion; however, iMaeda et al.l (120101 ) emphasize 
the importance of the viewing angle on the observed properties of Type la SNe. We cannot 
state with certainty whether SN 2001ay was a single or double-degenerate explosion. 

The A4130 line gives a velocity gradient of ?) = 226 ± 32 km s~^ d~^, while the A6355 
line gives v = 171 ± 35 km s~^ d~^. The Si II velo city gradient is considerably higher than 
70 km s~^ d~^, the criterion of iBenetti et al.l (120051 ) to include SN 2001ay with other "high 
velocity gradient" Type la SNe. 



Discussion 



4.1. Reddening 



As mentioned above, the Galactic reddening along the line of sight to SN 2001ay is 
E{B — V) = 0.026 mag, and the host-galaxy component is E{B — V) = 0.072 mag. Thus, for 
SN 2001ay, E{B — \^)totai ~ 0.098 m ag. Normal Galactic dust is charac terized by an average 
value of Rv = Ay / E{B — V^) = 3.1 (ICardelli. Clayton, fc Mathislll989l ). bu t dust associated 



with Type la SNe is often characterized by a lower value; see, for example, iKrisciunas et al. 
(120071 ). If we adopt Ry = 3.1 for the Galactic component of reddening for SN 2001ay and 
Rv = 2.4 ± 0.2 for the host-galaxy reddening, it follows that Ab = 0.35 ± 0.08, Ay = 
0.25 ± 0.06, Ar = 0.21 ±0. 05, Aj = 0.15 ± 0.04, Aj = 0.08 ± 0.02, and Ah = 0.04 ± 0.01 



mag (IKrisciunas et al.ll2006l Table 8 



4.2. Absolute Magnitudes at Maximum Light 



Garnavich et al.l ( 120041 ) give BVI decline-rate relations for w hat was then the known 



range of Ami5(_B) for Type la SNe. The slowest decliner used by iPrieto. Rest, fc Suntzefl 



^^However, iPoznanski et al.l (|201lh recently found that the uncertainty of the reddening derived from the 
Na I D Hues may be substantiahy larger than this. 
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( l2006bl ) for their light-curve fitting template algorithm was SN 1999aa, with Ami;^(B)= 0.83 
mag. For such a decline-rate parameter, the implied absolute magnitudes at maximum are 
Mb = -19.42, Mv = -19.39, and M/ = -18.85 on an Hq = 72 km s"^ Mpc"^ distance 
scale. Three of th e four super-Chandrasekhar-mass Type la SN candidates discussed by 



Scalzo et al.l (120101 ) are significantly more luminous than this. 



In Figure [To] we show the absolute \^-band magnitudes at maximum light of SNe 2001ay, 
2003fg, 2005eq, 2 006gz , 2007if, and 2009dc, along with the l^-band decline-rate relation of 



Garnavich et al.l (120041 ). Figure [TTl gives the near-IR absolute magnitudes of Type la SNe at 
their respective maxima, including SNe 200 lay and 2005eq. Three of the four super-Chandra 
candidates are overluminous, while all other very slow decliners have "normal" maximum 
brightness. 

Based on our best optical photometry (using PSF magnitudes and image-subtraction 
templates), we find apparent magnitudes for SN 2001ay of -Bmax = 16.71, \4iax = 16.64, 
-Rmax = 16.65, and /max = 16.79. These values include the K-corrections and S-corrections. 
Given the extinction values listed above and the distance modulus of the host galaxy given in 
Table m at maximum light we find Mb = -19.19 ± 0.12, My = -19.17, Mr = -19.10, and 
Ml = —18.90 mag (with uncertainties of ±0.10 mag); thus, SN 2001ay was noi overluminous 
at optical wavelengths. 

Typically, Type la S Ne peak in the near-IR ~ 3 d prior to the time of 5-band maximum 



(iKrisciunas et al.ll2004bl ). SN 2005eq peaked at least 7 d prior to B maximum. We surmise 
that SN 2001ay may have been ~ 0.12 mag brighter at maximum light than our first Js 
measurement, or Js(max) ^ 16.85 mag. In the H band SN 2001ay exhibited a very flat light 
curve; we adopt i/max = 16.97 ± 0.08 mag. The corresponding extinction-corrected absolute 
magnitudes at IR maximum are Mj^ = —18.77 ± 0.10 and Mfj = —18.62 ± 0.10. These 
values may be compare d to the mean values of all but the late-peaking fast decliners from 



Krisciunas et al.l (120091 ). namely (Mj) = —18.61 and (Mh) = —18.31 mag. The standard 
deviations of the distributions of the near-IR absolute magnitudes are about ±0.15 mag. 
In the Js and H bands, SN 200 lay was overluminous by ~ la and 2a, respectively; hence, 
SN 2001ay was not statistically significantly brighter in the near-IR than other Type la SNe. 



4.3. Bolometric Light Curve and Mass Budget 

In Figure [TTl we show the bolometric light curve of SN 2001ay, based on our broad-band 
photometry. Though Type la SNe near maximum light have spectral energy distributions 
that peak at optical wavelengths, they also emit UV and IR light. One might scale the 
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integrated flux by a factor of ~ l-l. a ccount for UV a nd IR light not i nclude d in the 
optical bandpasses (IStritzinger et al.ll2006l ). but according to lBlinnikov et al.l (|2006[ ) this still 
underestimates the bolometric flux. Adopting distance modulus 35.55 mag {D ^ 129 Mpc; 
see Tabled]) and a scale factor of 1.15, we obtain a peak derived bolometric luminosity {AttD"^ 
times the bolometric flux) of 1.20 x 10^^ erg s~^. Figure [T7I also shows the bolometric light 
curves of SNe 2007if and 2009dc. 

At maximum light the luminosity produced by radioactive ^^Ni is given by 



Ni 



where Ej^a is the energy input from the decay of ^^Ni, evalu ated at the time of bolometric 
maximum. Arnett's Rule implies that a = 1 (jArnettlll982l ). meaning that the gamma-ray 
deposition matches the bolometric flux at maximum light. However, the value of a can 
actually range from 0.8 to 1.3, depe nding on the explosion model. Fo r a delayed detona- 
tion model a = 1.2 is app r opriate (jKhokhlov. Miiller. fc HoflichI Il993l . Fig. 36); see also 
Hoflich fc Khokhlovl Jl996[ ). Iflowell et al.l tood ) adopt a = 1.2. 



In Figure ITS] we show the radioactive decay energy deposition function fit to the last few 
points of the bolometric light curve. In this figure we also show the cases of complete trapping 
of the 7 rays and complete 7-ray escape. We adopt a rise t ime of 18 d from ex plosion to bolo- 
metric maximum, cornparab le to observational results of iGarg et al.l (120071 ) for other Type 
la SNe. iHayden et al.l ( I2OIOI ) find an average rise time in t he B band of 17.38 d, with a range 
of 13 to 23 d, while the S-band rise time determined by iGaneshalingam. Li. fc Filippenko 
feoilh is about 18 d. Figure [18] uses a = l.oH 

The energy deposited by 1 Mq of ^^Ni into the explosion of a Type la SN is given by 
Stritzinger &: LeibundgutI ( 120051 ) as 



Em{l Mq) = (6.45 X 10^^)e"*«/'-^ + (1.45 x 10^^)e 



tR/8.8 



.43 N 



(2) 



where tn is the rise time in days from the moment of explosion to the time of bolometric 
maximum. The e-folding times of ^^Ni and ^^Co are 8.8 and 111.3 d, respectively. For tn = 



18 d and a = 1.0, Ln 



2.07 xl0^3 erg s"^ M"^; for a 



1.2, 



2.48 xl0^3 erg s"^ 



^■^The writers of this section could not co me to a consensus o n the best value of a to adopt, so in what 
follows we give the results for a — 1.0 and 1.2. 1 Baron et all (|2011l l choose a — 0.9 for their SN 2001ay model. 
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There was a nondetection of SN 2001ay on April 5.40 which was 17.1 rest-frame days 
before the time of 5-band maximum. The upper hmit of the brightness on that date was 
R ^ 19.5 mag, or ~ 3 mag fainter than the observed i?-band maximum. Although we have 
an image of the host galaxy at t = — 17 d, we would need to reach a much fainter magnitude 
limit to place a useful constraint on the time of the explosion. The lack of premaximum 
photometry does not allow us to determine the rise time of SN 2001ay. 

To obtain the ^^Ni yield in solar masses, we simply divide the peak luminosity (1.20 x 



43 



X a 



]^q43 gj^g g-i^ ii^Q coefficient given above for the adopted rise time (e.g., 2.07 x 10 
for in = 18 d); the result is 0.58/a Mq of ^^Ni. SN 2001ay certainly did not produce more 
than a Chandrasekhar mass of ^®Ni. 

A 5% uncertainty in the value of Hq leads to a 10% uncertainty in the luminosity 
calculated from the optical photometry. Given the uncertainties in the parameter a in 
Equation 1 and uncertainty in the adopted bolometric rise time, the minimum uncertainty 
in the ^^Ni yield is about ±0.15 M© divided by a. 



We would like to estimate the mass in the ejecta of SN 2001ay. Following I Jeffrey! (119991 ) 
and lStritzinger et al.l (120061 ) . the "fiducial time" (to) is the time scale for the ejecta to become 
optically thin: 



57r 



1 

2 1 



(3) 



where Mej is the ejecta mass, k is the 7-ray mean opacity, g is a dimensioi iless scale factor, 
and V e is the e-folding velocity of an exponential model's density profile. IStritzinger et al. 
(I2OO6I ) adopt K = 0.025 cm^ g~^ and Ve = 3000 km s~^. They also choose q = 1/3, meaning 
that ^^Ni was distributed throughout the ejecta; q = 1 for high concentrations of ^®Ni at the 
center of the ejecta. Using the last three points in our bolometric light curve (from t = 49 to 
111 d), adopting a = 1.2, and using a ^^Ni yield of 0.48 M0 we derive a "fiducial time" of 66.2 
± 3.0 d. Taking a = 1.0 and a ^^Ni yield of 0.58 Mq gives to = 57.2 ± 2.4 d. Both of these 



values are considerably longer than the 22 d to 35 d in Table 1 of IStritzinger et al 



Scalzo et al 



(2006), 
J2010f ). 



but not much greater than the 51 d for SN 2007if given in Figure 9 of 
The SN 2001ay fiducial time depends critically on the final point in the bolometric light 
curve, a point derived from PSF photometry without template subtractior i (hence possibly 
erroneous). Using to = 66.2 d and the values of q, Ve, and k adopted by IStritzinger et al. 



(I2OO6I ). the implied ejecta mass is 4.4 Mq, which is clearly wrong. For to = 57.2 d and the 



24 



Swift fc Li 



(|200l[) incorrectly give the date of this nondetection as 2001 April 4.4. 
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higher ^^Ni yield, the implied ejecta mass is 3.3 Mq. Adopting q = 1 instead could cut these 
down by a factor of three. Choosing g = 1/3 but = 1500 km s~^ would give 0.8 to 1.1 Mq 
for the ejecta mass. The lower e-foldi ng velocity is only slightly l ess than what one derives 
from the three-dimensional models of iRopke &: Hillebrandtl (120041 ). Given the uncertainties 
of the parameters necessary for Eq. 3 and the lack of high-quality photometry from ~ 50 to 
100 d after maximum light, we feel that a robust calculation of the ejecta mass of SN 200 lay 
is beyond the scope of this paper. 



We note that iTaubenberger et al.l (120111 ) derived an ejecta mass for SN 2009dc of 
2.8 M0, based on the expansion velocity and the timescale around maximum brightness, 
which depends on the optical opacity. They describe their result as "an utmost challenge for 
all scenarios that invoke thermonuclear explosions of white dwarfs." 



4.4. Spectroscopic Comparison and Modeling 



Wagers. Wang. &: AsztalosI (120 lOl ) have used wavelet analysis to search for spectroscopic 
correlations of Type la SNe. Their results are based on the analysis of a few dozen relatively 
nearby objects {z < 0.04). SN 2001ay is anomalous in a number of their plots, particularly 
their Figure 12, showing the correlation of the strength of the 4570 A emission feature vs. 
Ami5(5). 



In Figure [19] we illustrate a comparison of the spectra of SNe 200 lay, 2005eq, and 
2009dc. The most obvious difference between SN 2001ay and these other two objects is 
the much larger blueshift of Si II seen in SN 2001ay prior to B maximum. Also, the C II 
absorption in SN 2009dc is much stronger. SN 2005eq was apparently the hottest of the 
three, given the stronger presence of doubly ionized iron and t he weakness of sin gly ionized 
species. Thi s is more like the classi cal slow decliner SN 1991T ( iFilippenko 1119971 ). which is, 
in fact, how IContreras et al.l (I2OIOI Table 1) classify SN 2005eq. 



In Table [10] we give a summary of some observational characteristics of SNe 200 lay, 
2005eq, and 2009dc. Given the divergence of the light curves of SN 2005eq seen in Figure [6] 
compared to the ot her two objects, t he lar ger decline-rate parameter of SN 2005eq (Amis (5) 



0.78 mag) fro m IContreras et a. 



(I2OIOI) is more likely correct than the value of Ami5(S) 



0.72 mag from iFolatelli et al.l ( l2010l ). though we have used the latter to plot SN 2005eq 



in Figure [TUl SNe 2001ay and 2009dc have almost the same optical decline rate. SN 2001ay 
had normal optical peak brightness and high Si II velocity. SN 2009dc was overluminous at 
optical maximum and had much lower Si II velocity. Their H-hand light curves are unlike 
those of any other Type la SNe observed thus far, except for SN 1999aw. 



-17- 



To in vestigate the detai ls of our SN spectra, we use the SN spectrum-synthesis code 
SYNOW (IFisher et aLl 119971 ) . Although SYNOW has a simple, parametric approach to 
creating synthetic spectra, it is a powerful tool to aid line identifications which in turn provide 
insights into the spectral formation of the objects. To generate a synthetic spectrum one 
inputs a blackbody temperature (Tbb), a photospheric velocity (fph), and for each involved 
ion, an optical depth at a reference line, an excitation temperature (Texc), and the maximum 
velocity of the opacity distribution (fmax)- Moreover, it assumes that the optical depth 
declines exponentially for velocities above fph with an e-folding scale of v^. The strengths 
of the lines for each ion are determined by oscillator strengths, and the approximation of a 
Boltzmann distribution of the lower level populations is set by the temperature Tgxc- 

In Figure [2Ul we present our —1 d spectrum of SN 2001ay with a synthetic spectrum 
generated from SYNOW. This fit has Tbb = 20,500 K and Wph = 12,000 km s'^ The 
majority of the observed features are well fit by the synthetic spectrum. The ions used 
in the fit, as lab eled in the figure, are commonly observed in the near-maximum spectra 
of Type la SNe (IBranch et al.ll2005l ) with the exception of C II. Although the inclusion of 
C II in the fit produces an absorption feature at ~ 6300 A (due to C II A6580) that is too 
strong compared to the observed spectrum, it nonetheless yields a better fit to the absorption 
feature at ~ 4500 A (due to C II A4745), so we conclude that the inclusion of C II marginally 
improves the SYNOW fit. 

We attempt SYNOW fits to the SN 2001ay spectra at t ^ +6 d and 23 d. Figure E] 
shows the result for the t ^ +6 d spectrum. The model spectrum has Tbb = 14,000 K and 
fph = 11,000 km s~^, and has all the regular ions observed in a Type la SN. C II is no longer 
needed, but a relatively strong Na I line is now observed at ~ 5600 A. A SYNOW fit for the 
t 23 d spectrum (not shown) requires Tbb = 8500 K and fph = 9000 km s~^ 



Fi nally, in Figure 



I we show a model fit to our Keck spectrum using the SYNAPPS 
code of iThomas et al.l (|201l[) . The physical assumptions SYNAPPS uses match those of 



SYNOW ( IFisher et al.l 119971 ). so findings are restricted to identification of features and not 
quantitative abundances. But where SYNOW is completely interactive, SYNAPPS is auto- 
mated. This relieves the user from iterative adjustment of a large number of parameters (over 
50 variables) to gain fit agreement and assures more exhaustive searching of the parameter 
space. SYNAPPS can be thought of as the hybridization of a SYNOW-like calculation with a 
parallel optimization framework, where spectral fit quality serves as the objective function to 
optimize. A good fit constrains explosion models through interpretive spectral feature iden- 
tification, with the main result being the detection or exclusion of specific chemical elements. 
The velocity distribution of detected species within the ejecta can also be constrained. 



SYNAPPS indicates that SN 2001ay had a photospheric velocity of 10,800 km s ^ at 
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t = —1 d, which should be robust, as it is derived from 11 ions, namely C II, O I, Mg II, 
Si II, Si III, S II, Ca II, Fe II, Fe III, Co II, and Ni II. We are confident of the presence of 
C II. 



5. Conclusions 



We have presented the available spectra, as well as optical and near-IR photometry, of 
SN 2001ay. 

While SN 2001ay is the most slowly declining Type la SN ever discovered, it was not 
overluminous in optical bandpasses. In near-IR bands it was overluminous only at the 2a 
level at most. Unlike other very slow decliners such as SNe 2003fg, 2007if, and 2009dc, 
which were significantly overluminous, we do not have any evidence that SN 2001ay was a 
super- Chandrasekhar-mass explosion. 

SN 2001ay showed evidence for C II, but it was much weaker than in SN 2009dc. At 
early times Mg II and Si II were observed in SN 2001ay at high velocity (14,000 km s~^ 
and higher). By contrast, spectra of SN 2009dc did not show large velocities for Mg II and 
Si II. On the basis of a small number of super-Chandrasekhar-mass candidates, it seems that 
these objects exhibit rather low velocities, possibly the result of retardation due to a shell of 
material arising from the disruption of the less massive white dwarf in a double-degenerate 
system. 

SN 2001ay produced (0.58 ± 0.15)/a M© of ^^Ni, considerably less than a Chan- 
drasekhar mass. The value of a probably lies in the range 1.0 to 1.2. Naively, one might 
conclude that SN 2001ay was a single-degenerate explosion, but this is hardly a firm con- 
clusion. The very broad light curve might be explained by the trapping of the 7 rays in the 
explosi on. An explanatio n of the extremely slow decline will be discussed in a subsequent 



paper (IBaron et al.ll201ll ). 
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Table 1. Properties of SN 2001ay and its Host Galaxy 



J. tXl dliiCuC'l 


\/ Q 1 1 1 p 
V dl LLC 


Host galaxy 


IC 4423 


Host-galaxy type'^ 


Sbc 


Heliocentric radial velocity^ 


9067 km s"^ 


CMB-frame radial velocity^ 


9266 km s"^ 


Distance modulus*^ 


35.55 ± 0.1 mag 


E{B - V)Gal 


0.026 ± 0.006 mag 


E{B - V)host 


0.072 ± 0.008 mag 


SN a(J2000) 


14'^26'^17'!0 


SN (5(J2000) 


+26°14'55'.'8 


Offset from nucleus 


10'.'3W 9'.'3N 


Julian Date of i?-band maximum 


2,452,022.49 ±0.8 


UT Date of S-band maximum 


2001 April 23.0 


Ami5(E) 


0.68 ± 0.05 mag 


Mb (max) 


-19.19 ±0.12 


My (max) 


-19.17 ±0.10 


Mr (max) 


-19.10 ±0.10 


M/(max) 


-18.90 ±0.10 


Mj (max) 


-18.77 ±0.10 


Mh (max) 


-18.62 ±0.10 



Trom HyperLEDA. 



TOeers et al.l f ll995h . via NED. 
^From NED. 



'^Using Ho = 72 km s"^ Mpc^^ flFreedman et al.ll200l[ ) 
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Table 2. Optical Field-Star Sequence near SN 2001ay 



ID^ 


U (mag) 


B (mag) 


V (mag) 


R (ma 


g) 


/ (mag) 


1 


17.112 


(0.043) 


16.399 


(0.019) 


15.480 


(0.008) 


14.859 


(0 


018) 


14.386 


(0.015) 


2 


18.224 


(0.047) 


17.262 


(0.020) 


16.193 


(0.007) 


15.467 


(0 


019) 


14.921 


(0.016) 


3 


16.705 


(0.022) 


16.846 


(0.019) 


16.367 


(0.008) 


15.981 


(0 


022) 


15.622 


(0.021) 


4 


17.933 


(0.035) 


18.271 


(0.030) 


17.864 


(0.005) 


17.503 


(0 


025) 


17.131 


(0.025) 


5 


19.246 


(0.174) 


18.729 


(0.027) 


17.884 


(0.015) 


17.302 


(0 


024) 


16.796 


(0.033) 


6 


16.813 


(0.031) 


16.853 


(0.017) 


16.320 


(0.007) 


15.917 


(0 


017) 


15.574 


(0.014) 


7 


20.437 


(0.129) 


19.585 


(0.015) 


18.313 


(0.010) 


17.487 


(0 


017) 


16.730 


(0.010) 


8 


20.137 


(0.230) 


19.987 


(0.196) 


18.674 


(0.070) 


18.110 


(0 


100) 


17.467 


(0.111) 


9 


21.270 


(0.289) 


19.802 


(0.108) 


18.790 


(0.059) 


18.009 


(0 


101) 


17.467 


(0.074) 


10 














20.719 


(0 


155) 


18.930 


(0.061) 



*^Magnitude uncertainties (la) are given in parentheses. 

''The IDs are the same as in Fig. [TJ Star 10, located 18.4" SE of Star 7, is not visible 
in the l^-band finder. 



Table 3. Infrared Field Star Sequence near SN 2001ay 



ID^ 


Js (mag) 


H (mag) 


1 


13.718 (0.014) 


13.251 (0.011) 


6 


15.153 (0.021) 


14.881 (0.008) 


7 


15.869 (0.007) 


15.281 (0.007) 


8 


16.566 (0.030) 


15.932 (0.031) 


10 


17.838 (0.024) 


17.359 (0.029) 



^Magnitude uncertainties (Icr) are 
given in parentheses. 

^The IDs are the same as in Fig. [1] 
and Table [21 
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Table 4. Fully Corrected Optical Photometry of SN 2001ay^ 





U (mag) 


B (mag) 


V (mag) 


R (mag) 


I (mag) 


Telescope"^ 


Subs'* 


2020.78 




16.773 


^0.043) 


16.708 (0.010) 


16.694 (0.015) 


16.771 


(0.017) 


1 


N 


2021.71 


16.419 (0.034) 


16.731 


(0.020) 


16.700 (0.010) 


16.695 (0.006) 


16.802 


(0.010) 


1 


N 


2022.90 




16.687 


(0.050) 


16.643 (0.020) 


16.635 (0.040) 


16.783 


(0.050) 


2 


Y 


2023.88 




16.724 


(0.070) 


16.640 (0.020) 


16.661 (0.030) 


16.781 


(0.020) 


2 


Y 


2024.89 




16.709 


(0.040) 


16.647 (0.020) 


16.626 (0.020) 


16.829 


(0.030) 


2 


Y 


2025.89 




16.737 


(0.030) 


16.666 (0.020) 


16.692 (0.020) 


16.868 


(0.050) 


2 


Y 


2026.86 




10. * DD 


^{J.KJjUJ ) 


16.664 (0.020) 


16.659 (0.030) 


lO.OOO 


[U.U4U j 


2 


Y 


2028.89 




16.841 


(0.100) 


16.699 (0.080) 








2 


Y 


2030.88 




16.913 


(0.040) 


16.734 (0.020) 


16.775 (0.030) 


16.937 


(0.040) 


2 


Y 


2032.88 




17.065 


(0.050) 


16.799 (0.030) 


16.860 (0.030) 


16.952 


(0.050) 


2 


Y 


2037.88 




17.290 


1^0.050) 


16.939 (0.030) 


16.948 (0.020) 


17.102 


(0.040) 


2 


Y 


2039.86 




1 7 c;q c; 
i / .OOO 


^U.UoU ) 


16.945 (0.050) 


16.963 (0.040) 


1 ( .yjo ( 


((\ n'7n\ 
[U.U ( U ) 


2 


Y 


2041.80 




17.D47 


J 0.040) 


17.089 (0.050) 








2 


Y 


2043.85 








17.143 (0.020) 


16.968 (0.020) 


17.093 


(Ci ^^^^ 


2 


Y 


2045.86 




17.952 


(0.050) 


17.266 (0.030) 


17.006 (0.020) 


17.024 


(0.040) 


2 


Y 


2047.84 




18.134 


(0.040) 


17.298 (0.040) 


17.015 (0.030) 


16.969 


(0.040) 


2 


Y 


2051.81 




18.250 


1^0.050) 


nAU (0.030) 


17.123 (0.020) 


16.902 


(0.040) 


2 


Y 


2059.81 




18.554 


(0.080) 


17.794 (0.040) 


17.385 (0.020) 


17.138 


(0.040) 


2 


Y 


2064.83 




18.773 


(0.080) 


18.064 (0.070) 


17.629 (0.040) 


17.251 


(0.060) 


2 


Y 


2069.79 




1 Q not; 


^v.yj (yj) 


18.117 (0.060) 


17.813 (0.030) 


i / .ouo 


(r\ nftn\ 

i^U.UDU ) 


2 


Y 


2074.74 




1 Q n7i 


^U.U t U) 


18.219 (0.070) 


18.008 (0.030) 


± / .oUO 


(ft f\af\\ 

i^U.UDU ) 


2 


Y 


21)75.81 




1Q QQQ 


L U.UUU J 


18.153 (U.UbUj 


18.099 (U.U5Uj 


17 604 




2 


Y 










1 c ocrn /n ^'7^^ 
lo.z5U (U.UrUj 


18.18U (U.UoUJ 




(r\ ^^^^^ 


2 


Y 






Id. ( Uy 


( n nQ'3^ 
^U.Uoo ) 


Id. 655 (U.013) 


16.693 (0.011) 


Id. / oo 


in m ftA 
(^U.UlD ) 


3 


Y 










16.629 (0.025) 


16.686 (0.011) 


1 fi 7R1 




3 


Y 


2029.81 




16.958 


(0.060) 


16.716 (0.010) 


16.781 (O.OlOj 


16.910 


(0.015) 


3 


Y 


2UOD.C5D 




17.322 


(0.020) 


1 OQO ni n\ 
iD.yoz (U.UlUj 


1 c\ nni /n ^OQ^ 
lo.yUl yj.yJZo) 


17.007 


(0.024) 


Q 
O 


IN 


2037.83 


17.503 (0.040) 


17.348 


(0.019) 


16.990 (0.010) 


16.967 (0.018) 


17.099 


(0.020) 


3 


N 


2038.77 




17.488 


(0.009) 




17.043 (0.010) 


17.079 


(0.010) 


3 


N 


2040.79 


17.655 (0.031) 


17.578 


(0.014) 


17.054 (0.010) 


16.964 (0.016) 


16.939 


(0.016) 


3 


N 


2041.84 




17.668 


(0.034) 


17.003 (0.020) 


16.949 (0.033) 


16.910 


(0.034) 


3 


N 


2043.78 




17.784 


(0.010) 




16.990 (0.010) 


16.845 


(0.017) 


3 


N 


2046.82 










17.107 (0.010) 


16.871 


(0.010) 


3 


N 


2050.82 




18.269 


(0.013) 




16.992 (0.085) 


16.816 


(0.010) 


3 


N 


2063.76 








18.083 (0.017) 




17.225 


(0.029) 


3 


N 


2140.70 




20.110 


(0.130) 


19.520 (0.100) 


19.390 (0.080) 


19.030 


(0.090) 


3 


N 


2026.74 




16.841 


(0.062) 


16.673 (0.040) 








4 


N 


2030.67 




16.953 


(0.017) 


16.725 (0.010) 








5 


N 


2046.62 




18.063 


(0.024) 


17.194 (0.013) 








5 


N 


2047.65 




18.209 


(0.047) 


17.247 (0.027) 








5 


N 


2055.63 








17.687 (0.028) 








5 


N 


2052.65 


18.633 (0.041) 


18.356 (0.023) 


17.468 (0.012) 


17.152 (0.027) 


16.850 (0.029) 


6 


Y 



''The BVRI data in this tabic include the K-corrections, S-corrcctions, and corrections applied to photometry derived 
using aperture magnitudes without subtraction templates. Magnitude uncertainties (Icr) are given in parentheses. 

''Julian Date minus 2,450,000. 
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ci = CTIO 0.9-m, 2 = KAIT, 3 = Nickel 1-m, 4 = Yale-AURA-Lisbon-Ohio (YALO) 1-m, 5 = LCO 2.5-m, and 6 

= CTIO 1.5-m telescopes. 

''image-subtraction templates used? Y = yes; N = No. 
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Table 5. K-Corrections for SN 2001ay in Optical Bands'^ 



UT Date (2001) 


Phase (d)'' 


Telescope 


AS (mag) 


AV (mag) 


AiJ (mag) 


A/ (mag) 


Apr 21 


-2 


FLWO 1.5-m 


-0.010 


-0.031 


-0.120 


-0.045 


Apr 24 


+ 1 


FLWO 1.5-m 


0.000 


-0.028 


-0.118 


-0.068 


Apr 26 


+3 


MMT 6.5-m 


-0.001 


-0.020 


-0.090 


-0.029 


Apr 27 


+4 


FLWO 1.5-m 


0.018 


-0.014 


-0.097 


-0.015 


Apr 30 


+7 


FLWO 1.5-m 


0.010 


-0.017 


-0.094 


-0.050 


May 02 


+9 


FLWO 1.5-m 


0.023 


-0.005 


-0.093 


-0.065 


May 10 


+17 


BAG 2.16-m 




-0.024: 


-0.064 




May 16 


+23 


FLWO 1.5-m 


0.080 


0.047 


-0.055 


0.042 


May 23 


+30 


FLWO 1.5-m 


0.090 


0.066 


-0.057 


0.030 


May 24 


+31 


MMT 6.5-m 


0.097 


0.077 


-0.039 


-0.008 


May 25 


+32 


FLWO 1.5-m 


0.085 


0.068 


-0.061 


0.025 


May 25 


+32 


MMT 6.5-m 


0.100 


0.081 


-0.040 


0.022 


May 30 


+37 


FLWO 1.5-m 


0.082 


0.070 


-0.049 


0.019 


Jun 18 


+56 


FLWO 1.5-m 


0.089 


0.073 


-0.067 


0.043 


Jun 18 


+56 


MMT 6.5-m 


0.093 


0.070 


-0.056 


0.059 



^These corrections (measured in magnitudes) are subtracted f rom the photo metric data to correct the 
photometry to the rest frame. These were calculated using the iBesselll lll990l ) filter profiles. 

''The number of observer-frame days since the time of B maximum, 2001 April 23. 
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Table 6. Interpolated Corrections to Optical Photometry' 







Kcv 




Kr, 


B 




R 




Telescope^ 


2020.78 


-0.010 


-0.031 


-0.120 


-0.053 


-0.008 


-0.003 


-0.006 


-0.002 


1 


2021.71 


-0.006 


-0.030 


-0.116 


-0.052 


-0.007 


-0.004 


-0.005 


-0.006 


1 


2022.90 


-0.005 


-0.028 


-0.113 


-0.050 


-0.008 


-0.005 


0.022 


-0.017 


2 


2023.88 


-0.002 


-0.025 


-0.110 


-0.048 


-0.008 


-0.005 


0.021 


-0.017 


2 


2024.89 


0.004 


-0.022 


-0.106 


-0.046 


-0.007 


-0.005 


0.020 


-0.017 


2 


2025.89 


0.006 


-0.021 


-0.103 


-0.044 


-0.007 


-0.005 


0.019 


-0.016 


2 


2026.86 


0.007 


-0.019 


-0.100 


-0.041 


-0.007 


-0.005 


0.019 


-0.016 


2 


2028.89 


0.013 


-0.013 






-0.006 


-0.004 






2 


2030.88 


0.022 


-0.006 


-0.088 


-0.033 


-0.005 


-0.002 


0.017 


-0.016 


2 


2032.88 


0.030 


0.000 


-0.083 


-0.028 


-0.005 


-0.001 


0.017 


-0.016 


2 


2037.88 


0.045 


0.016 


-0.071 


-0.017 


-0.005 


0.005 


0.017 


-0.015 


2 


2039.86 


0.050 


0.022 


-0.066 


-0.012 


-0.005 


0.007 


0.017 


-0.015 


2 


2041.80 


0.057 


0.029 






-0.006 


0.008 






2 


2043.85 




0.037 


-0.059 


-0.007 




0.010 


0.019 


-0.014 


2 


2045.86 


0.071 


0.045 


-0.056 


0.002 


-0.007 


0.011 


0.020 


-0.014 


2 


2047.84 


0.078 


0.053 


-0.054 


0.007 


-0.008 


0.011 


0.021 


-0.014 


2 


2051.81 


0.089 


0.068 


-0.050 


0.015 


-0.011 


0.012 


0.023 


-0.013 


2 


2059.81 


0.089 


0.075 


-0.048 


0.031 


-0.017 


0.009 


0.027 


-0.011 


2 


2064.83 


0.089 


0.073 


-0.050 


0.039 


-0.018 


0.007 


0.029 


-0.010 


2 


2069.79 


0.087 


0.067 


-0.053 


0.045 


-0.018 


0.004 


0.030 


-0.009 


2 


2074.74 


0.088 


0.066 


-0.058 


0.048 


-0.011 


0.005 


0.030 


-0.007 


2 


2075.81 


0.088 


0.067 


-0.059 


0.049 


-0.009 


0.005 


0.030 


-0.007 


2 


2079.75 




0.072 


-0.062 


0.050 




0.010 


0.028 


-0.006 


2 


2022.99 


-0.005 


-0.028 


-0.113 


-0.050 


0.003 


-0.009 


0.024 


0.027 


3 


2024.91 


0.004 


-0.022 


-0.106 


-0.046 


0.004 


-0.008 


0.023 


0.036 


3 


2029.81 


0.017 


-0.010 


-0.091 


-0.035 


0.000 


-0.002 


0.022 


0.060 


3 


2036.86 


0.042 


0.013 


—0.073 


—0.019 


—0.014 


0.007 


0.024 


0.073 


3 


2037.83 


0.045 


0.016 


-0.071 


-0.017 


-0.016 


0.008 


0.024 


0.076 


3 


2038.77 


0.047 




-0.068 


-0.014 


-0.018 




0.024 


0.079 


3 


2040.79 


0.053 


0.026 


-0.064 


-0.010 


-0.022 


0.011 


0.024 


0.082 


3 


2041.84 


0.057 


0.029 


-0.063 


-0.007 


-0.024 


0.011 


0.024 


0.082 


3 


2043.78 


0.064 




-0.059 


-0.003 


-0.028 




0.024 


0.082 


3 


2046.82 






-0.056 


0.004 






0.025 


0.075 


3 


2050.82 


0.087 




-0.051 


0.013 


-0.037 




0.027 


0.053 


3 


2063.76 




0.074 




0.037 




-0.005 




0.000 


3 


2026.74 


0.007 


-0.019 














4 


2030.67 


0.008 


-0.006 














5 


2046.62 


0.075 


0.049 














5 


2047.65 


0.078 


0.053 














5 


2055.63 




0.075 














5 


2052.66 


0.089 


0.068 


-0.050 


0.017 


0.018 


0.006 


0.000 


- 0.028 


6 



''Values in Columns 2 through 9 are in magnitudes. The K-corrrections are subtracted from the data in 
Table 4. The S-corrections are added to the data in Table |4] 

Julian Date minus 2,450,000. 

= 1 = CTIO 0.9-m, 2 = KAIT, 3 = Nickel 1-m, 4 = Yale-AURA-Lisbon-Ohio (YALO) 1-m, 5 = LCD 
2.5-m, and 6 = CTIO 1.5-m telescopes. 
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Table 7. Near-Infrared Photometry of SN 2001ay^ 



JD'^ Js (mag) H (mag) Telescope^ 



2024 


73 










17 


045 


(0 


042) 


1 


2025 


67 


16 


977 


(0 


032) 


16 


893 


(0 


044) 


1 


2026 


70 


17 


050 


(0 


028) 










1 


2029 


70 


17 


132 


(0 


022) 


17 


000 


(0 


035) 


1 


2031 


68 


17 


136 


(0 


022) 


16 


965 


(0 


036) 


1 


2032 


70 


17 


190 


(0 


028) 


17 


110 


(0 


027) 


1 


2045 


66 


17 


291 


(0 


024) 


16 


873 


(0 


020) 


2 


2047 


65 


17 


509 


(0 


049) 










1 


2049 


66 


17 


396 


(0 


035) 


16 


853 


(0 


039) 


1 


2051 


65 


17 


430 


(0 


034) 


16 


894 


(0 


038) 


1 


2053 


67 


17 


337 


(0 


031) 










1 


2055 


67 


17 


387 


(0 


031) 


17 


078 


(0 


056) 


1 


2056 


67 


17 


420 


(0 


031) 


16 


998 


(0 


022) 


1 


2061 


65 










16 


905 


(0 


022) 


2 


2068 


57 


17 


592 


(0 


033) 


17 


157 


(0 


028) 


2 



^Magnitude uncertainties (Icr) are given in parenthe- 
ses. 

'^Julian Date minus 2,450,000. 

^1 = LCO 1-m, 2 - LCO 2.5-m telescopes. 



Table 8. Spectroscopic Observations in the "CfA Set" 



UT Date 


Julian Day 


Tel. 


Range 


Res. 


Pos.A. 


Par.A. 


Airmass 


Flux Std. 


Sky 


Seeing 


Slit 


Exp. 


Observer(s) 


(2001) 






(A) 


(A) 


n 


n 








(") 


(") 


(s) 




Apr 21.48 


2452020.98 


FLWO 


3720-7540 


7.0 


71.00 


70.02 


1.48 


Feige34 




1-2 


3 


1200 


Calkins 


Apr 24.31 


2452023.81 


FLWO 


3720-7540 


7.0 


0.00 


-16.25 


1.01 


Feige34 


clear 


2 


3 


1200 


Berlind 


Apr 26.33 


2452025.83 


FLWO 


3720-7540 


7.0 


45.00 


43.51 


1.01 


Felge34 


some clouds 


1-2 


3 


1200 


Berlind 


Apr 26.40 


2452025.90 


MMTO 


3250-8900 


8.0 


Par.A. 


71.15 


1.12 


Feige34/BD26 








900 


Matheson 


Apr 26.41 


2452025.91 


MMTO 


3250- 8900 


8.0 


Par.A. 


71.52 


1.15 


Feige34/BD26 








900 


Matlieson 


Apr 27.36 


2452026.86 


FLWO 


3720-7540 


7.0 


64.00 


65.79 


1.04 


HZ44 




1-2 


3 


1200 


Calkins 


Apr 30.36 


2452029.86 


FLWO 


3720-7540 


7.0 


70.00 


67.86 


1.05 


Feige34 


clear 


1-2 


3 


1200 


Berlind 


May 02.29 


2452031.79 


FLWO 


3720-7540 


7.0 


-9.00 


-19.75 


1.01 


HZ44 


clear 


1-2 


3 


1200 


Berlind 


May 16.26 


2452045.76 


FLWO 


3720-7540 


7.0 


—33.00 


—13.29 


1.01 


HZ44 


clear 


1-2 


3 


1200 


Calkins 


May 23.33 


2452052.83 


FLWO 


3720-7540 


7.0 


69.00 


71.42 


1.14 


Feige34 


cirrus 


2 


3 


1200 


Rines 


May 24.32 


2452053.82 


MMTO 


3220-8900 


8.0 


at Par.A. 


70.92 


1.11 


BD28/BD26 


clear 






1200 


Challis, Phelps 


May 24.33 


2452053.83 


MMTO 


4950-10000 


8.0 


Par.A. 


71.51 


1.16 


BD26 


clear 






1200 


Challis, Phelps 


May 25.32 


2452054.82 


MMTO 


3220-8900 


8.0 


Par.A. 


71.36 


1.14 


Fcigc34/BD26 


clear 






1800 


Challis 


May 25.35 


2452054.85 


MMTO 


4950-10000 


8.0 


Par.A. 


71.52 


1.23 


BD26 


clear 






1800 


Challis ' 
Rmes InO 


May 25.37 


2452054.87 


FLWO 


3720-7500 


7.0 


69.00 


70.63 


1.39 


Feige34 


cirrus 


2 


3 


1200 


May 30.36 


2452059.86 


FLWO 


3720-7540 


7.0 


71.00 


70.51 


1.40 


HZ44 


clear 


1-2 


3 


1200 


Calkins 1 


May 30.38 


2452059.88 


FLWO 


3720-7540 


7.0 


71.00 


69.68 


1.53 


HZ44 


cirrus 


1-2 


3 


1200 


Calkins 


.luiic 18.22 


2452078.72 


FLWO 


3720-7540 


7.0 


69.00 


66.99 


1.05 


HZ44 


clouds 


1-2 


3 


1200 


Berlind 


June 18.29 


2452078.79 


MMTO 


3200-8900 


8.0 


Par.A. 


71.24 


1.29 


BD28/BD26 


some clouds 






900 


Matheson 


June 18.30 


2452078.80 


MMTO 


3200-8900 


8.0 


Par.A. 


70.84 


1.35 


BD28/BD26 


some clouds 






900 


Matheson 
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Table 9. Other Spectroscopic Observations 



UT Date (2001) 


Phase (d)'' 


Telescope 


Wavelength Range (A) 


Resolution (A) 


Apr 22 


-1 


Keck II 10-m 


3922-10000 


1.25 


Apr 29 


+6 


KPNO 2.1-m 


5350-[9500] 


4.8 


Apr 30 


+7 


Lick 3-m 


3260-10600 


5.5/10.5 


May 01 


+8 


LCO 2.5-m 


3500-9199 


7.0 


May 02 


+9 


HST 2.4-m 


2877-5348 


3.55 to 4.10 


May 09 


+16 


HST 2.4-m 


2930-5695 


3.55 to 4.10 


May 10 


+17 


BAO 2.16-m 


4100-9000 


9.6 


May 16 


+23 


Lick 3-m 


3300-10450 


5.5/10.5 


May 17 


+24 


LCO 2.5-m 


3500-9350 


7.0 


May 20 


+27 


BAO 2.16-m 


3450-9000 


9.6 


May 23 


+30 


KPNO 4.0-m 


3500-5493/5300-11600 


7.8/13.8 


May 25 


+32 


BAO 2.16-m 


3700-9300 


9.6 


May 26 


+33 


LCO 2.5-m 


3800-9234 


7.0 


May 30 


+37 


BAO 2.16-m 


3970-8200 


4.8 



''Number of observer frame days since time of S-band maximum, 2001 April 23. 



Table 10. Comparison of Three Objects 



Parameter 


SN 200 lay 


SN 2005eq 


SN 2009dc 


Ami5(5) 


0.68 (0.05) 


0.72-0.78 


0.72 (0.03) 


My (max) 


-19.17 (0.10) 


-19.40 (0.08) 


-20.16 (0.10) 


Si II velocity^ (km s"^) 


-14,340 [-1.5 d] 


-10,200 [-5 d] 


-8600 [-7 d] 


C II line(s) 


weak 


not detectable 


very strong 


if-band light curve 


flat 


peaked early 


increasing 



'^For the Aq = 6355 A line. The values in brackets correspond to the numbers 
of rest-frame days after the time of B maximum w hen the line was measu red 
at this velocity. For SN 2009dc the velocity is from [Silverman et al.l (120101 ). 
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Fig. 1. — Finder chart for IC 4423, SN 200 lay, and some field stars in our Galaxy. 

Fig. 2. — K-corrections for BVRI photometry of SN 2001ay. 

Fig. 3. — S-corrections for B- and y-band photometry of SN 2001ay. 



Fig. 4. — Effective J-band transmission functions for the iBesselll (119901 ) filter, the KAIT 
filter, and the Nickel 1-m filter. 

Fig. 5. — S-corrections for R- and /-band photometry of SN 2001ay. 

Fig. 6. — Comparis on of the B - and l^-band light curves of SNe 200 lay, 2005eq 
dContreras eTaPboiol ). and 2009dc ^Silverman et al.lboioh . 

Fig. 7. — Unreddened B — V colors of SN 200 lay vs. time since B-band maximum. 
Fig. 8. — Similar to Figure O except for the R and / bands. 



Fig. 9.— Near-IR pho tometry of SNe 2009dc dStritzinger et al.ll201lh . 200 lay, and 2005eq 
JContreras et al.lboioh . 

Fig. 10. — Absolute \^-band magnitudes at maximum light of six very slowl y declining Type 



la SNe vs. the decline-rate parameter Ami^{B). The decline-rate relation of lGarnavich et al. 
(I2OO4J ) is also shown. 

Fig. 11. — Near-IR absolute magnitudes at maximum light of Type la SNe. 
Fig. 12. — Temporal sequence of spectra of SN 200 lay. 

Fig. 13. — Additional spectra of SN 2001ay obtained with the Lick Observatory 3-m Shane 
telescope, HST, and the Las Campanas Observatory 2.5-m du Pont telescope. 

Fig. 14. — A portion of our highest signal-to-noise ratio spectrum of SN 2001ay, taken with 
the Keck ESI on 2001 April 22. 

Fig. 15. — Profile of the Na I D lines in our Keck ESI spectrum. The Milky Way and 
host-galaxy components are clearly present. 

Fig. 16. — Blueshifted velocity of two Si II lines in the spectra of SN 2001ay around maximum 
hght. 



Fig. 17.— Bolometric lum inosities (L) of SNe 2001ay, 2007if flScalzo et al.ll2010[ ). and 2009dc 
(jTaubenberger et"aL 2011 ). measured in erg s~^ 
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Fig. 18. — The best-fit radioactive decay energy deposition function (middle solid line) to 
the UVOIR light curve (blue squares) of SN 2001ay. 

Fig. 19. — Comparison of the optical spectra of SNe 2001ay, 2005eq, and 2009dc. 

Fig. 20. — SYNOW modeling of the optical spectrum of SN 2001ay at 1 d before B maximum. 
Upper curve = data; lower curve = SYNOW model. The dashed lines indicate the addition 
of C II to the model. 

Fig. 21. — SYNOW modehng of the optical spectrum of SN 2001ay 6 d after B maximum. 
Upper curve = data; lower curve = SYNOW model. 

Fig. 22. — SYNAPPS modeling of the optical spectrum of SN 2001ay 1 d before B maximum. 
The blue dashed line (model spectrum) can be shifted arbitrarily along the vertical axis to 
match the actual spectrum (shown in black). 
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Krisciunas et al. Figure [H Finder chart for IC 4423, SN 2001ay, and some field stars in our 
Galaxy. This 7.4' x 7.4' image is a 200 s \^-band exposure obtained on 2001 May 23 with 
the CTIO 1.5-m telescope. North is up and east is to the left. The SN is located 10.3" W 

and 9.3" N of the galaxy core. 
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Days since B-band maximum 

Krisciunas et al. Figure [2l K-corrections for BVRI photometry of SN 200 lay. These 
account for the shifting of the SN spectrum to longer wavelengths owing to the redshift of 
the host galaxy. These values are to be subtracted from the photometry. Values derived 
from the high signal-to-noise ratio MMT spectra are plotted as larger symbols. 
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Krisciunas et al. Figure [3l S-corrections for B- and V"-band photometry of SN 2001ay, 
based on spectra from the FLWO 1.5-m and the MMT. The individual points are shown 
only for the KAIT corrections. The corrections for the older Nickel 1-m chip are essentially 
the same as those shown for the newer Nickel 1-m chip. S- corrections are added to the 
photometry to correct it to the filter system of iBesselll (119901 ). 
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Krisciunas et al. Figure HI Effective J-band transmission functions for tlie iBesselll ( 1l990h 

filter, the KAIT filter, and the Nickel 1-m filter used with the newer, higher quantum 
efficiency chip. Each filter profile has been normalized to 1.00 at maximum throughput. 
We also include the t = +23 d spectrum from Figure [13] (multiplied by 10^^); it shows 
significant flux past 9000 A, which is measured by the Lick 1-m system but excluded in the 

KAIT measurements. 
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Krisciunas et al. Figure [51 S-corrections for R- and /-band photometry of SN 2001ay. 
Corrections for the older Nickel 1-m chip are shown as a dotted line; corrections for the 
new Nickel 1-m chip are shown as a dashed line. For the /-band corrections with the Nickel 
1-m filter and both CCD chips used in the camera, we can only use spectra that extend to 
sufficiently long wavelengths (our Keck spectrum, two spectra taken with the Lick 3-m, and 
one with the KPNO 4-m). For the Nickel 1-m photometry to match the KAIT photometry, 
we would need the S-corrections to be a factor of 3 larger than the deri ved valu es. 



S-corrections are added to the photometry to correct it to the filter system of iBesselll (1l990l ). 
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Krisciun as et al. Figure E Co mparison of th e B- and l^-band ligh t curves of SNe 200 lay, 



I 



2005eq fIContreras et al.l 120101 ). and 2009dc fISilverman et al.l 120101 ). The SN 2001ay data 
are K-corrected and S-corrected, and the CTIO 0.9-m and Nickel 1-m aperture photometry 
includes additional corrections discussed in the text. For SN 200 lay the upward pointing 
triangles represent data which result from PSF photometry using image-subtraction 
templates. The downward pointing triangles are SN 200 lay data based on aperture 

photometry. 
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Krisciunas et al. Figure [71 Unreddened B — V colors of SN 2001ay vs. time since S-band 
maximum. Upward pointing triangles are data derived using image-subtraction templates. 

Downward pointing triangles represent data derived fro m aperture magnitudes without 
template subtractions. The dashed line is the locus from lPrieto. Rest, fc Suntzefa (l2006bl ) 

for their most slowly declining object (with A mi^f - B) = 0.83 mag). The soli d line is the 
commonly used Lira-Phillips locus (jLiralll995l : iPhillips et al.lll999l ). 
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Krisciunas et al. Figure [H Similar to Figure |6l except for the R and I bands. In the case 
of SN 2005eq the data are taken in Sloan r and i filters, not in Kron-Cousins R and /. The 

Nickel 1-m J-band data are not shown. 
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Krisciunas et al. Figu re [9l Near-IR photome try of SNe 2009dc (IStritzinger et al.ll201ll ). 
2001ay, and 2005eq (IContreras et al.ll2010l ). The sohd hne is an educated guess that 
SN 2001ay may have been ~ 0.12 mag brighter in Js at maximum hght than our earhest Js 

datum. 




Krisciunas et al. Figure [TOl Absolute V^-band magnitudes at maximum light of six very 
slowly declining Type la SNe vs. the decline-ra t e par ameter Ami5(-B). The decline-rate 
relation of iGarnavich et al.l (120041 ) is also shown. 
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Krisciunas et al. Figure [TTl Near-IR absolute magnitudes at maximum light of Type la 
SNe. While SNe 2001ay and 2005eq are more luminous than the average of other objects in 
the H band, one would not consider them significantly overluminous. The H-hand point 
for SN 2009dc corresponds to the same epoch as in the J-band plot, but SN 2009dc 
brightened 0.26 mag in H over the following month. The diamond-shape points correspond 
to objects that peak in the near-I R later than B raaximu m. These objects are subluminous 
in all bands. See iKrisciunas et al.l (l2009l ) for more details. 
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Krisciunas et al. Figure [121 Temporal sequence of spectra of SN 2001ay. The labels 
indicate the number of observer-frame days with respect to 5-band maximum. Eleven 

which have no telescope indicated were taken with the FLWO 1.5-m. The t = +6 d 
spectrum from the KPNO 2.1-m telescope has been combined with the t = +9 d HST 
spectrum. The solid vertical line marks the nominal wavelength of the Aq = 6355 A line of 
Si II. The dashed vertical line shows Si II line blueshifted by 10,000 km s^^. Telluric 

oxygen is identified by an Earth symbol. 
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Krisciunas et al. Figure fT3l Additional spectra of SN 2001ay obtained with the Lick 
Observatory 3-m Shane telescope, HST, and the Las Campanas Observatory 2.5-m du Pont 
telescope. The captioning within the figure and the meaning of the vertical lines are the 
same as in Figure [T2j Telluric oxygen is identified by an Earth symbol. 
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Krisciunas et al. Figure [HI A portion of our highest signal-to-noise ratio spectrum of 
SN 2001ay, taken with the Keck ESI on 2001 April 22. The dashed lines for Si III and S II 
correspond to blueshifts of 9000 km s~^. The dashed magenta lines for C II correspond to a 
blueshift of 12,000 km s~^. The dot-dashed red lines for Mg II and Si II correspond to a 
blueshift of 14,000 km s~^. The dotted line for Na I is the rest wavelength. 
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Krisciunas et al. Figure [151 Profile of tlie Na I D lines in our Keck ESI spectrum. The 
Milky Way and host-galaxy components are clearly present. 
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Krisciunas et al. Figure [161 Blueshifted velocity of two Si II lines in the spectra of 
SN 200 lay around maximum light. Within la the gradient derived from the two lines is 
the same, about 200 km s~^ d"^. This qualifies SN 200 lay to be a "high velocity gradient" 

object jBenetti et aliboosh . 
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Krisciunas et al. Fi gure [T71 Bolometric lumin osities (L) of SNe 200 lay, 2007if (jScalzo et al. 
20101 ). and 2009dc ( iTaubenberger et al.ll201lh . measured in erg s~^. For SN 2001ay we have 
scaled our integrated bolometric luminosity by a factor of 1.15 to account for flux not 
included in our photometric bandpasses. We adopted a distance modulus of 

M = 35.55 mag from Tabled] 
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Krisciunas et al. Figure [181 The best-fit radioactive decay energy deposition function 
(middle solid line) to the U VOIR light curve (blue squares) of SN 200 lay. Here we adopt 
Arnett's Rule (jArnettlll982l ). with a = 1.0, which stipulates that the gamma-ray deposition 
matches the bolometric flux at maximum light. The top solid line corresponds to the case 
of complete trapping of 7 rays and positrons (i.e., r ^ 1), while the bottom solid line is the 

case of complete 7-ray escape (r ^ 1). The vertical dashed line indicates the epoch at 
which the fit begins. The best fit to the latest few points suggests that the ejecta became 
optically thin to 7 rays 57.2 ± 2.4 d past maximum brightness. 
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Krisciunas et al. Figure UHl Comparison of the optical spectra of SNe 2001ay, 2005eq, and 

2009dc. 
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Krisciunas et al. Figure l20l SYNOW modeling of the optical spectrum of SN 2001ay at 1 d 
before B maximum. Upper curve = data; lower curve = SYNOW model. The dashed lines 

indicate the addition of C II to the model. 
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Krisciunas et al. Figure EH SYNOW modeling of the optical spectrum of SN 2001ay 6 d 
after B maximum. Upper curve = data; lower curve = SYNOW model. 
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SN 2001ay (t = -ld) 
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Krisciunas et al. Figure |22l SYNAPPS modeling of the optical spectrum of SN 2001ay 1 d 
before B maximum. The blue dashed line (model spectrum) can be shifted arbitrarily 
along the vertical axis to match the actual spectrum (shown in black). 



